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We present initial results of a numerical investigation of the chiral symmetry restoring transition in the (3+1)- 
dimensional Nambu - Jona-Lasinio model with both non-zero baryon chemical potential (/is 7^ 0) and isospin 
chemical potential {fii 7^ 0). With non-zero isospin chemical potential, the model suffers from a sign problem. 
We proceed in two ways: (i) We perform "partially quenched" simulations in which /i/ is made non-zero only 
during the measurement of chiral observables; (ii) We perform full simulations with imaginary isospin chemical 
potential with the aim to analytically continue results to real /ij. 



1. Motivation 

At asymptotically high baryon chemical poten- 
tial (fJ-s) and low temperature (T) - where QCD 
can be treated in a perturbative manner - the 
ground-state of quark matter is found to be that 
of a colour-superconductor (for a recent review 
see e.g. ^). Unfortunately, the sign problem pre- 
vents us from using lattice QCD to determine the 
ground-state at the more moderate densities typ- 
ical in the cores of compact (neutron) stars. 

One way to proceed is to study model field the- 
ories such as the Nambu - Jona-Lasinio (NJL) 
model. This purely fermionic field theory, in 
which colour-neutral quarks interact via a four- 
point contact term, not only contains the same 
global symmetries as two flavour QCD, but can 
be simulated on the lattice even with /is 7^ 0. 
In [2] we show that the ground-state of the lat- 
tice model with = l^d = fJ-B, i-e. with "up" 
and "down" quarks sharing a common Fermi sur- 
face, exhibits s-wave superfluidity via a standard 
BCS pairing between quarks of different flavours; 
i.e. 



{ud) ^ 0; Ascs ^ 0. 



(1) 



Within the cores of compact stars, however, the 
Fermi momenta /c]i. and kp are expected to be 



separate. A simple argument based on that of |3j 
suggests that for a two flavour Fermi liquid of 
massless quarks and electrons with /is = dOOMeV 
and both weak equilibrium (/i^ = + /ig) and 
charge neutrality (2n„/3 — n^/S — = 0) en- 
forced, all the Fermi momenta of the system are 
determined: 
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^J'u = ^J'B — Me/2 = 355.5MeV, 
fid^tJ-B+ = 444.5MeV, 



(2) 



89MeV. 



The effect of separating the Fermi surfaces of 
pairing quarks in QCD should be to make the 
colour-superconducting phase less energetically 
favourable. Introducing fj,j oc (/x^ — fid) 7^ could 
prove a good method, therefore, to investigate the 
stability of the superfluid phase. 

2. The Model 

The action of the lattice NJL model (with a 
1) is given by 



INJL 



(3) 
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where ^ = (u, d)'"' is the SU{2) doublet of stag- 
gered up and down quarks defined on lattice sites 
X and ^ = a + iir.f is a matrix of bosonic aux- 
iliary fields defined on dual sites x. The fermion 
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kinetic matrix M,j.y is defined in O and we choose 
the same bare parameters used therein. 

One can separate the Fermi surface of up and 
down quarks by simultaneously setting baryon 
chemical potential = (Mji + ['^d)!'^ 7^ ^^^d 
isospin chemical potential /i/ = — Mrf)/2 7^ 0. 
With ^1 = 0^ T2MT2 = M*, which is a sufficient 
condition to show that det M is both real and pos- 
itive With /i/ 7^ however, this is no longer 
true such that once again we are faced with the 
sign problem. 

The fact that physically the two scales are or- 
dered iJLi suggests that one may be able 
to apply techniques recently developed to study 
QCD with iiB <- T |51CT7) . First, however, we 
present the results of a partially quenched calcu- 
lation. 

3. Partially Quenched /i/ 

Whilst the primary motivation for investigat- 
ing /i/ ^ is to study the superfluid phase which 
sets in at large /is , this requires one to introduce 
an explicit symmetry breaking parameter (j); it 
is currently not clear how to study /i/ 7^ in the 
j ^ limit 2 . Instead, we choose to study the 
chiral symmetry restoring phase transition with 
the aim of controlling the systematics of intro- 
ducing /i/ ^ 0. 

The first step we take is to perform a "partially 
quenched" calculation in which fij — when gen- 
erating the background fields and is made non- 
zero only during the measurement of fermion ob- 
servables. In particular, we measure the up and 
down quark condensates 

as functions of /is for various /i/ on a 12^ lattice. 
Some results are presented in Fig. 

The results agree qualitatively with those of 
mean-field studies of the model in which the in- 
troduction of a small but non-zero /z/ is seen to 
suppress the up quark condensate and enhance 
the down |8I9| . This can be understood by not- 
ing that with fij 7^ 0, /i„ (fid) reaches the critical 
chemical potential at a lower (higher) value oi ijb- 
Why (uu) deviates from the /.ij = result more 
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Figure 1. (uu) and ^dd) condensates for various 
Hb and 11 1 on a 12^ lattice. 

than (dd) is not so clear, but may be an effect 
associated with T > on a finite lattice 

4. Continuation from Imaginary pij 

One recent method used to study QCD with 
small baryon chemical potential fiB is to measure 
quantities at imaginary /i^, where the measure of 
the path integral is real, and fit results to a trun- 
cated Taylor expansion in /is/T about fiB =^ 0. 
One can then analytically continue the results to 
real ^b [Z|- We propose to use a similar method 
by simulating at imaginary isospin chemical po- 
tential (/i/), where once again detM is real, and 
expanding observables in powers of flj //is . 

For our initial investigation we have measured 
the quark condensates and baryon and isospin 
number densities 

1 d\nZ 1 - 

"^■^^279^^4^"^°"^'^^"'^'^ 

on a 12* lattice at hb = 0.6, which from Fig.^can 
be seen to be where the effect of having fij ^ 
is largest. In QCD, one can show that e.g. (i'lp) 
expanded about ^b = is analytic in such 
that for small imaginary ^b the quantity remains 
real For our simulations, however, this is not 
the case, and measured quantities are, in general, 
complex. Therefore, we fit the data by the Taylor 
series 

(|i:;j)^|(t)(0' 
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Figure 2. Real and imaginary parts of (uw) and 
(dd) as functions of imaginary fij with /is = 0.6 
on a 12^ lattice. 

and 

(-)^|(S;)(f^)^ 

each truncated at some suitable point. We can 
then analytically continue to real /i/ using e.g. 

{uu) = Ao + Aim - Ml^] - Asfi^j + ■■■ (8) 

Figure |2 shows the real and imaginary parts 
of the condensates as functions of jl] with fits to 
constant + quadratic and linear only forms re- 
spectively. The quality of the fits is very good, 
suggesting that our ansatze are correct. The real 
parts are in complete agreement whilst the imagi- 
nary parts are anti-correlated, which implies that 
the sum of the condensates {iptp} contains only 
the even terms of the expansion, whilst their dif- 
ference {^T^ip') contains only the odd terms. A 
similar effect is seen in the number densities, as 
riB is found to be real and therefore even in 
M//M-Bi whilst the isospin density is pure imag- 
inary. For a good quality fit to nj, however, we 
have to include the cubic term. These data and 
fits are plotted in Fig. O 
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Figure 3. Baryon and isospin density as functions 
of imaginary fij with jiB = 0.6 on a 12-^ lattice. 



Whilst these results are only preliminary, we 
have shown that we can calculate the coefficients 
in © and (0 as functions oi jiB and in principle, 
reproduce reliable forms of the curves in Fig. ^ 
With this aim, we plan to repeat this exercise 
for various values of /is in both the the chirally 
broken and restored phases on various lattice vol- 
umes. Also, whilst it is difficult to study the di- 
quark condensate in the j — > limit, it would be 
interesting to compare the response of {ud) to /x/ 
at fixed j to that of l^ipi') at fixed mass. 
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